The Astrophysical Journal Letters, 756:L29 (4pp), 2012 September 10 
© 2012. The American Astronomical Society. All rights reserved. Printed in the U.S.A. 


doi: 1 0. 1 088/204 1 -8205/756/2/L29 


SIZE DISTRIBUTIONS OF SOLAR FLARES AND SOLAR ENERGETIC PARTICLE EVENTS 

E. W. Oliver 1 , A. G. Ling 2 , A. Belov 3 , and S. Yashiro 4,5 

1 Space Vehicles Directorate, Air Force Research Laboratory, Sunspot, NM 88349, USA 
3 Atmospheric Environmental Research, Lexington, MA 02421, USA 
3 IZMIRAN, Troitsk, Moscow Region 142190, Russia 
4 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 
5 Department of Physics, The Catholic University of America, Washington, DC 20064, USA 
Received 2012 July 10; accepted 2012 August 7; published 2012 August 21 

ABSTRACT 

We suggest that the flatter size distribution of solar energetic proton (SEP) events relative to that of flare soft X-ray 
(SXR) events is primarily due to the fact that SEP flares are an energetic subset of all flares. Flares associated with 
gradual SEP events are characteristically accompanied by fast (ft: 1 000 km s _1 ) coronal mass ejections (CMEs) 
that drive coronal/interplanetary shock waves. For the 1996-2005 interval, the slopes (a values) of power-law size 
distributions of the peak 1-8 A fluxes of SXR flares associated with (a) >10 MeV SEP events (with peak fluxes 
^1 pr cm~ 2 s _1 sr _1 ) and (b) fast CMEs were ~ 1.3-1. 4 compared to ~ 1.2 for the peak proton fluxes of > 10 MeV 
SEP events and ~2 for the peak 1-8 A fluxes of all SXR flares. The difference of ~0. 15 between the slopes of the 
distributions of SEP events and SEP SXR flares is consistent with the observed variation of SEP event peak flux 
with SXR peak flux. 
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1. INTRODUCTION 

Hudson (1978) pointed out that the size distribution of 
solar energetic proton (SEP) events is flatter than those of 
flare electromagnetic emissions. When size distributions are 
represented by power laws, of the form 

dN/dx ~ x - “ events/unitx, (1) 

SEP events are found to have values of a that range from ~1.2 
to 1.4 versus values of ~1.8 to 2.0 for flare electromagnetic 
emissions. Table 1 contains representative reported values of 
a for the two types of emissions. The SEP distributions are 
relatively independent of the energy range considered and the 
flare source longitude (Belov et al. 2007). 

The cause of the difference between the slopes of the dis- 
tributions in Table 1 is not well understood. Schrijver et al. 
(2012) recently listed four possible contributors to the discrep- 
ancy, including “(1) SEP spectral distributions may depend on 
event energy (which could include a dependence on the parti- 
tioning between flare radiative and coronal mass ejection, CME, 
bulk-kinetic energies), (2) background corrections, (3) effects of 
compound events involving two or more CME/shocks . . ., and 
(4) particle propagation effects in the heliosphere.” Schrijver 
et al. (2012) considered a fifth, geometrical, possibility, the di- 
lution of SEP fluence over an opening angle into the heliosphere 
and calculated that this effect might be primarily responsible for 
the observed drop off or “kink” in the >10 MeV SEP fluence 
size distribution at a level of ~5 — 10 x 10 9 pr cm" 2 first reported 
by Lingenfelter & Hudson (1980). 

As an alternative explanation for the difference in the power- 
law slopes (values of a) of the distributions in Table 1, Hudson 
(1978) suggested that peak proton flux ( J P ) varied with the total 
flare energy (W; assumed to scale with peak soft X-ray, SXR, 
flux) as J P ~ WP with ft > 1 . A determination of ft from the 
known values of a in Table 1 for SXR events (a x = 1.84; Drake 
1971) and SEP events (a P = 1.15 ± 0.05; Van Hollebeke et al. 


1975) using the formula 

ft = («x — 1)/(«p - 1) (2) 

yielded ft = 5.6(+2.8, —1.4). Hudson suggested that “proton 
flares probably do not represent a separate category of so- 
lar flares, but stand out principally because proton emission 
strongly biases the observer toward the more energetic events.” 

In the same year as Hudson’s paper, Kahler et al. (1978) pre- 
sented the first direct evidence that SEP events were associated 
with CMEs, indicating that proton flares did, in fact, represent a 
separate class of flares which in time came to be known as erup- 
tive flares to distinguish them from compact or confined events 
that lacked CMEs (Svestka & Cliver 1992). Current thinking is 
that large “gradual” SEP events are caused by fast CMEs that 
drive coronal/interplanetary shocks manifested by radio type 
II bursts (Reames 1999; Gopalswamy et al. 2002; Cliver et al. 
2004). In this paper, we use this modern view of SEP flares 
to revisit the cause of the difference in the slopes of the size 
distributions of solar flare electromagnetic emissions and SEP 
events. 

Our analysis is presented in Section 2 and results are summa- 
rized and discussed in Section 3. 

2. ANALYSIS 

From lists of SEP events for solar cycle 23 compiled by 
Belov et al. (2005, 2007) and Cane et al. (2010), we made 
a list of all >10 MeV proton events from 1996 to 2005 that 
were associated with flares with nominally good (W20-W85) 
magnetic connection to Earth. To focus on gradual SEP events, 
we limited the sample to events with hourly averaged peak 
fluxes ^1 pr cm -2 s~ sr _1 (or A 1 proton flux unit (pfu)). 
Cliver (2009) suggested that the upper limit peak flux for the 
smaller “impulsive” SEP events is ~3 pfu at this energy (see 
Cliver & Cane 2002 for a discussion of SEP event terminology). 
We only considered the prompt component of each SEP event, 
ignoring delayed peaks associated with interplanetary shocks. 
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Figure 1. Size distributions for (1) peak 1-8 A fluxes of all ^Ml.O SXR flares 
(black line); (2) peak 1-8 A fluxes of SXR flares associated with >10 MeV 
proton events with peak fluxes ^1 pfu (light blue line); (3) peak 1-8 A fluxes 
of SXR flares associated with ^1000 km s _1 CMEs (red line); and (4) peak 
proton fluxes of >10 MeV SEP events (reference distribution for slope only, 
magenta line). The SXR flare classification is given at the top of the figure. All 
distributions are based on flares or SEP events that originated in a longitude 
range from W20 to W85 (1996-2005). 


Table 1 

Values of a for SEP and Flare Size Distributions 


Parameter 

Value of a 

Reference 

SEPs 

20-80 MeV Flux 

1.15 ± 0.05 

Van Hollebeke et al. (1975) 

>10 MeV Fluence 

1.24 + 0.04 

Gabriel & Feynman (1996) 

>10 MeV Flux 

1.37 + 0.03 

Belov et al. (2005) 

Flare radiation 
Hard X-rays 

1.73 ± 0.01 

Crosby et al. (1993) 

Microwaves 

1.73 ± 0.04 

Nita et al. (2004) 

Soft X-rays 

1.98 + 0.11 

Aschwanden & Freeland (2012) 


We subtracted the pre-event background for both the 1-8 A SXR 
bursts and >10 MeV SEP events (extrapolating the time profile 
of any preceding SEP event to the time of the peak of the event). 
The list of 58 SEP events with their associated SXR and CME 
data is given in Table 2. 

Figure 1 contains power-law size distributions for (1) the 
peak 1-8 A fluxes of 540 SXR flares of class ^Ml.O (^l.Ox 
10 -5 W m -2 ) (black line; data from Yashiro et al. 2006, 
http://cdaw.gsfc.nasa.gov/pub/yashiro/flare_cme/fclist_pub.txt); 
(2) the peak 1-8 A fluxes of 52 ^ M1.0 SXR flares associated 
with >10 MeV proton events (light blue line; Table 2); (3) the 
peak 1-8 A fluxes of 59 SXR flares associated with CMEs with 
speeds Js 1 000 km s _1 (red line; Yashiro et al. 2006); and (4) the 
peak proton fluxes of 58 >10 MeV SEP events (reference distri- 
bution for slope only, magenta line; Table 2). The 1000 km s' 1 
CME speed threshold corresponds roughly to the speeds of 
CMEs required to produce type II bursts that will extend from 
the metric to the kilometric wavelength range (Gopalswamy 
et al. 2005). Such shocks from western-hemisphere flares are 
highly associated with large SEP events (Gopalswamy et al. 
2008). 

To facilitate comparisons, and minimize SEP propagation 
and SXR occultation effects, all of the size distributions in 
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Figure 2. Scatter plot of peak > 10 MeV SEP event flux vs. flare peak 1-8 A flux 
for proton events originating from W20 to W85 heliolongitude (1996-2005), 
with geometric mean regression line. 


Figure 1 are based on flares or SEP events that originated 
in a longitude range from W20 to W85 for the period from 
1996 to 2005. Distributions 2-4 are based on background (pre- 
event) corrected SXR and SEP data. No background subtraction 
was made for distribution 1 but the a value of 2.10 ± 0.12 
is comparable to that (1.98 ± 0.11) recently obtained in 
the comprehensive study of Aschwanden & Freeland (2012; 
Table 1). The smaller slope (1.73 ± 0.01) of the hard X-ray flare 
distribution (Crosby et al. 1993) is attributed to the fact that big 
flares produce relatively more hard X-rays (e.g., Battaglia et al. 
2005). The slopes of the distributions based on uncorrected data 
(for size distributions 2 (1.28 ± 0.24), 3 (1.38 ± 0.22), and 
4 (1.15 ± 0.12)) differ little from those of the background- 
corrected distributions shown in Figure 1 . 

Figure 1 shows that the slopes of size distributions 2 and 
3, for the peak 1-8 A fluxes of SEP- and fast-CME-associated 
flares, respectively, are significantly flatter than that for the peak 
1-8 A fluxes of all flares. The similarity of distributions 2 and 3 
reflects substantial overlap in the two databases; 61% (36/59) of 
the fast-CME-flares were associated with SEP events in Table 1 
and 75% (41/55) of the SEP flares (no CME coverage for the 
three events) in Table 1 had associated fast CMEs. The median 
CME speed of the 14 events with speeds < 1000 km s 1 was 
806 km s -1 (range from 464 to 998 km s -1 ). The a values of 
these two distributions (SEP flares: a — 1.31 ± 0.24; fast CMEs: 
1.39 ± 0.22) are close to that of the > 10 MeV SEPs during this 
period (distribution 4: 1.16 ± 0.12). 

Inserting the value of ax = 1.31 ±0.24 from distribution 2 for 
the SXR peak fluxes of SEP flares into Equation (2) and using 
oip = 1 . 16± 0. 12 gives a value for /3 of 1.94 (+11.9, —1.69). The 
1 .94 base value of ft is similar to that obtained for the regression 
line in the scatter plot in Figure 2 between the peak proton 
fluxes of the 58 > 10 MeV SEP events in Table 2 and their 
associated peak 1-8 A SXR fluxes (1.59 ± 0.25; background 
subtracted for both parameters). The regression line in Figure 2 
was determined by assuming that the uncertainties in the logs 
of the >10 MeV flux and the SXR intensity are comparable, an 
assumption we believe is justified because the approximately 
three-order-of-magnitude scatter in CME energy for a given 
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Table 2 

SEP Events with Peak Fluxes ^1 pfu (1996-2005) 


Date 

(yyyy/mm/dd) 



SXR 



>10 MeV 

CME 

Onset 

Time 

(UT) 

Peak 

Time 

(UT) 

Initial 

Class 

BKG Corn 
Class 

Solar 

Lat/Lon 

Peak 

Flux 

(pfu) 

BKG Corr. 
Flux 
(pfu) 

Speed 
(km s -1 ) 

1997/11/04 

05:52 

05:58 

X 2.1 

X 2.1 

S14W33 

6.59e+01 

6.56e+01 

785 

1997/11/06 

11:49 

11:55 

X 9.4 

X 9.4 

S18W63 

4.69e+02 

4.68e+02 

1556 

1998/05/06 

07:58 

08:09 

X 2.7 

X 2.6 

S11W65 

1.52e+02 

1.52e+02 

1099 

1998/09/30 

13:08 

13:50 

M 2.8 

M 2.7 

N19W85 

9.79e+02 

9.79e+02 

d-g- 

1998/11/07 

11:02 

11:06 

M 2.4 

M 2.3 

N14W43 

6.46e+00 

5.76e+00 

632 

1998/11/22 

06:30 

06:42 

X 3.7 

X 3.7 

S27W82 

2.85e+00 

2.60e+00 

d-g- 

1999/06/04 

06:52 

07:03 

M 3.9 

M 3.8 

N17W69 

5.08e+01 

4.88e+01 

2230 

2000/02/12 

03:51 

04:10 

M 1.7 

M 1.6 

N26W24 

1.74e+00 

1.59e+00 

1107 

2000/04/04 

15:12 

15:41 

C 9.7 

C 8.1 

N16W66 

5.07e+01 

5.05e+01 

1188 

2000/05/15 

15:46 

16:01 

C 7.8 

C 4.7 

S22W68 

1.30e+00 

1.14e+00 

1212 

2000/06/10 

16:40 

17:02 

M 5.2 

M 5.0 

N22W40 

4.07e+01 

4.04e+01 

1108 

2000/06/18 

01:52 

01:59 

X 1.0 

M 9.9 

N23W85 

3.07e+00 

2.82e+00 

629 

2000/06/23 

14:18 

14:31 

M 3.0 

M 2.9 

N23W72 

1.31e+00 

1.15e+00 

847 

2000/06/25 

07:17 

07:52 

M 1.9 

M 1.8 

N16W55 

3.42e+00 

3.27e+00 

1617 

2000/07/22 

11:17 

11:34 

M 3.7 

M 3.3 

N14W56 

1.54e+01 

1.36e+01 

1230 

2000/11/08 

22:42 

23:28 

M 7.4 

M 7.0 

N10W75 

1.20e+04 

1.20e+04 

1732 

2001/01/28 

15:40 

16:00 

M 1.5 

M 1.3 

S04W59 

4.33e+01 

4.31e+01 

916 

2001/04/02 

10:58 

11:36 

X 1.1 

M 9.9 

N16W62 

3.17e+00 

2.40e+00 

992 

2001/04/02 

21:32 

21:51 

X20.0 

X19.9 

N17W78 

6.59e+02 

6.57e+02 

2505 

2001/04/12 

09:39 

10:28 

X 2.0 

X 2.0 

S20W42 

4.40e+01 

3.90e+01 

1184 

2001/04/15 

13:19 

13:50 

X14.4 

X14.4 

S20W84 

9.00e+02 

8.99e+02 

1199 

2001/04/26 

11:26 

13:12 

M 7.8 

M 7.7 

N17W27 

3.63e+00 

3.44e+00 

1006 

2001/06/04 

16:11 

16:33 

C 3.2 

C2.7 

N23W60 

1.91e+00 

1.58e+00 

464 

2001/09/15 

11:04 

11:28 

M 1.5 

M 1.3 

S24W51 

9.20e+00 

9.03e+00 

478 

2001/10/19 

16:13 

16:30 

X 1.6 

X 1.6 

N15W30 

9.68e+00 

5.28e+00 

901 

2001/11/22 

20:18 

20:36 

M 3.8 

M 3.7 

S24W68 

2.00e+01 

1.97e+01 

1443 

2001/11/22 

22:32 

23:30 

M 9.9 

M 9.5 

S15W34 

4.34e+03 

4.32e+03 

1437 

2001/12/26 

04:32 

05:40 

M 7.1 

M 6.6 

N08W54 

7.24e+02 

7.24e+02 

1446 

2002/02/20 

05:52 

06:12 

M 5.1 

M 4.9 

N12W72 

8.44e+00 

8.27e+00 

952 

2002/04/17 

07:46 

08:24 

M 2.6 

M 2.4 

S14W36 

1.76e+01 

1.75e+01 

1240 

2002/04/21 

00:43 

01:51 

X 1.5 

X 1.5 

S14W84 

2.32e+03 

2.32e+03 

2393 

2002/05/22 

03:18 

03:54 

C 5.0 

C 2.5 

S22W53 

1.08e+02 

1.08e+02 

1557 

2002/08/14 

01:47 

02:12 

M 2.3 

M 2. 1 

N09W54 

2.19e+01 

2.17e+01 

1309 

2002/08/20 

08:22 

08:26 

M 3.4 

M 3.2 

S11W38 

1.77e+00 

1.46e+00 

1099 

2002/08/22 

01:47 

01:57 

M 5.4 

M 5. 1 

S07W62 

3.13e+01 

3.12e+01 

998 

2002/08/24 

00:49 

01:12 

X 3.1 

X 3.1 

S02W81 

3.03e+02 

3.03e+02 

1913 

2002/11/09 

13:08 

13:23 

M 4.6 

M 4.3 

S12W29 

3.16e+02 

3.16e+02 

1838 

2003/05/28 

00:17 

00:27 

X 3.6 

X 3.6 

S06W20 

1.13e+01 

1.12e+01 

1366 

2003/05/31 

02:13 

02:24 

M 9.3 

M 9.2 

S07W65 

1.58e+01 

1.56e+01 

1835 

2003/10/26 

17:21 

18:19 

X 1.2 

X 1.2 

N02W38 

3.76e+02 

3.76e+02 

1537 

2003/11/02 

17:03 

17:25 

X 8.3 

X 8.2 

S14W56 

1.31e+03 

1.31e+03 

2598 

2003/11/03 

01:09 

01:30 

X 2.7 

X 2.6 

N10W83 

1.42e+03 

1.02e+03 

827 

2003/11/04 

19:29 

19:44 

X28.0 

X28.0 

S19W83 

3.03e+02 

3.02e+02 

2657 

2004/04/11 

03:54 

04:19 

C 9.6 

C 9.5 

S14W47 

2.70e+01 

2.68e+01 

1645 

2004/07/25 

14:19 

15:14 

M 1.1 

M 1.0 

N08W33 

6.76e+01 

6.74e+01 

1333 

2004/09/19 

16:46 

17:12 

M 1.9 

M 1.9 

N03W58 

4.58e+01 

4.56e+01 

d-g- 

2004/10/30 

16:18 

16:33 

M 5.9 

M 5.8 

N15W20 

1.94e+00 

1.04e+00 

690 

2004/11/09 

16:59 

17:19 

M 8.9 

M 8.8 

N07W51 

6.80e+01 

4.30e+01 

2000 

2004/11/10 

01:59 

02:13 

X 2.5 

X 2.5 

N09W49 

2.66e+02 

2.54e+02 

3387 

2005/01/17 

06:59 

09:52 

X 3.8 

X 3.8 

N15W25 

4.41e+03 

4.30e+03 

2094 

2005/01/19 

08:03 

08:22 

X 1.3 

X 1.1 

N19W47 

1.10e+02 

7.97e+01 

2020 

2005/01/20 

06:36 

07:01 

X 7.1 

X 7. 1 

N14W61 

1.44e+03 

1.41e+03 

882 

2005/05/06 

03:05 

03:14 

C 9.3 

C 8.7 

S06W74 

1.30e+00 

1.08e+00 

1120 

2005/05/06 

11:11 

11:28 

M 1.3 

M 1.3 

S05W80 

1.85e+00 

1.35e+00 

1144 

2005/07/09 

21:47 

22:06 

M 2.8 

M 2.7 

N11W27 

2.39e+00 

2.22e+00 

1540 

2005/07/13 

14:01 

14:49 

M 5.0 

M 4.9 

NI0W80 

1.02e+01 

9.99e+00 

1423 

2005/08/22 

00:44 

01:33 

M 2.6 

M 2.6 

S09W48 

5.63e+00 

5.41e+00 

1194 

2005/08/22 

16:46 

17:27 

M 5.6 

M 5.6 

S12W60 

3.11e+02 

3.1 le+02 

2378 
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peak SXR flux (Yashiro & Gopalswamy 2008) is similar to the 
range in the scatter of SEP peak flux with CME speed (Kahler 
2001). Applying a solar longitude correction factor (Kahler 
1982) to the peak SEP fluxes in Figure 2 had a negligible effect 
on the slope of the regression line; it yielded a fi of 1.61 ± 0.23. 

3. SUMMARY AND DISCUSSION 

We have presented evidence that the difference in the slopes 
of the power-law size distributions of solar flares and SEP 
events arises primarily because SEP flares, in particular those 
associated with gradual SEP events, represent an energetic 
subset of all flares, characterized by fast (Js 1 000 km s _1 ) CMEs. 
The slope (a value) of the size distribution for the peak 1-8 A 
fluxes of SEP- (or fast CME-) associated SXR flares is ~ 1.3—1 .4, 
compared to ~1.2 for the peak proton fluxes of >10 MeV 
SEP events and ~2 for the peak 1-8 A fluxes of all SXR flares 
(Figure 1). The remaining difference in the slopes of the peak 
flux distributions of SEP flares and SEP events is consistent 
with the observed variation of SEP peak flux with SXR peak 
flux (Figure 2). This second aspect of the analysis is poorly 
constrained, however, and all or part of the remaining difference 
could result from a propagation effect such as the dilution of SEP 
flux with opening angle proposed by Schrijver et al. (2012). In 
either case, it is clear that much of the discrepancy in the slopes 
of the size distributions of SEP events and SEP flares can be 
accounted for by the energetic nature of gradual SEP flares in 
relation to all flares. 

Yashiro et al. (2006) determined power-law slopes for dis- 
tributions of the peak 1-8 A fluxes of different classes of SXR 
flares as follows: flares without CMEs (a = 2.52 ± 0.03), all 
flares (2.16 ± 0.03), and flares with CMEs of any speed (1.98 ± 
0.05). The further decrease in a that we obtained for flares with 
Js 1 000 km s' 1 CMEs was expected. For type II bursts, Pearson 
et al. (1989) constructed histograms of the peak counting rates 
of hard X-ray bursts with and without associated metric type II 
bursts. These histograms have a ~ 1.2 for the type-II-associated 
events and ~2 for the non-type-II events. Gopalswamy et al. 
(2005) showed that a hierarchal relationship exists between the 
frequency range of type II bursts and the speeds/energies of as- 
sociated CMEs: (1) metric only type IIs have a mean CME 
speed of 610 km s -1 , (2) metric to decametric/hectometric 
type II bursts are associated with CMEs with a mean speed 
of 1115 km s -1 , and (3) metric to kilometric type IIs have a 
mean CME speed of 1490 km s -1 . Belov et al. (2007) were the 
first to draw attention to the similarity of the size distributions 
of SEP flares and SEP events. 

In regard to SEPs, there are three kinds of flares: (1) those that 
produce large “gradual” events (Reames 1999; Cliver 2009); 

(2) those that produce small “impulsive” SEP events with 
abundance anomalies for 3 He and high-Z elements (Reames 
et al. 1985; Mason et al. 1986; Reames & Ng 2004); and 

(3) confined flares that do not produce escaping particles. The 
inferred magnetic topologies of these three kinds of flares are 
quite different, involving the pinching off of a rising loop for 
gradual events (e.g., Forbes 2000), interchange reconnection 
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between open and closed field lines for impulsive events (e.g., 
Shimojo & Shibata 2000; Kahler et al. 2001), and field relaxation 
within a closed flux tube or interaction between closed flux tubes 
for confined events. What is remarkable is that these disparate 
types of flares combine to form a single well-defined power- 
law size distribution over several orders of magnitude in SXR 
intensity. 
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